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This is an introduction to models of accretion discs around black holes. After a presentation of 
the non-relativistic equations describing the structure and evolution of geometrically thin accretion 
discs we discuss their steady-state solutions and compare them to observation. Next we describe 
in detail the thermal-viscous disc instability model and its application to dwarf novae for which 
it was designed and its X-ray irradiated-disc version which explains the soft X-ray transients, 
i.e. outbursting black-hole low-mass X-ray binaries. We then turn to the role of advection in 
accretion flows onto black holes illustrating its action and importance with a toy model describing 
both ADAFs and slim discs. We conclude with a presentation of the general-relativistic formalism 
describing accretion discs in the Kerr space-time. 


I. INTRODUCTION 


The author of this chapter is old enough to remember the days when even serious astronomers doubted the existence 
of accretion discs and scientists snorted with contempt at the suggestion that there might be such things as black 
holes; the very possibility of their existence was rejected, and the idea of black holes was dismissed as a fancy of 
eccentric theorists. Today, some 50 years later, there is no doubt about the existence of accretion discs and black 
holes; both have been observed and shown to be ubiquitous in the Universe. The spectacular ALMA image of the 
protostellar disc in HL Tau [35] is breathtaking and we can soon expect to see the silhouette of a supermassive black 
hole in the center of the Galaxy in near infrared m or millimeter-waves GDI- 

Understanding accretion discs around black holes is interesting in itself because of the fascinating and complex 
physics involved but is also fundamental for understanding the coupled evolution of galaxies and their nuclear black 
holes, i.e. fundamental for the understanding the growth of structures in the Universe. The chance that inflows onto 
black holes are strictly radial, as assumed in many models, are slim. 

The aim of the present chapter is to introduce the reader to models of accretion discs around black holes. Because 
of the smallness of black holes the sizes of their accretion discs span several orders of magnitude: from close to the 
horizon up 100 000 or even 1 000 000 black-hole radii. This implies, for example, that the temperature in a disc 
around a stellar-mass black hole varies from 10' K, near the its surface, to ~ 10 3 K near the disc’s outer edge at 10 s 
black-hole radii, say. Thus studying black hole accretion discs allow the study of physical regimes relevant also in a 
different context and inversely, the knowledge of accretion disc physics in other systems such as e.g. protostellar discs 
or cataclysmic variable stars, is useful or even necessary for understanding the discs around black holes. 

Section [II] contains a short discussion of the disc driving mechanisms and introduces the a-prescription used in this 
chapter. In Section m after presenting the general framework of the geometrically thin disc model we discuss the 
properties of stationary solutions and the Shakura-Sunyaev solution in particular. The dwarf-nova disc instability 


model and its application to black-hole transient sources is the subject of Section IV The role of advection in accretion 
onto black holes is presented in Section [V] with the main stress put on high accretion rate flows. Finally, Section [VT| 
and |VII| about the general-relativistic version of the accretion disc equations concludes the present chapter. 
Notations and definitions 

The Schwarzschild radius (radius of a non-rotating black hole) is 


2 GM - M 

R s = —= 2.95 x 10 5 —— cm, 

CT M(7) 


(1) 


where M is the mass of the gravitating body and c the speed of light. 


* Chapter 1 in Astrophysics of Black Holes - From fundamental aspects to latest developments , Ed. Cosimo Bambi, Springer: Astrophysics 
and Space Science Library 440, (2016); DOI: 10.1007/978-3-662-52859-4. 
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The Eddington accretion rate is defined as 


Afedd 


^Edd 

rjC 2 


1 47 tGM 

V CK es 


1 2ncRs 

V K es 


1.6 x 10 18 r?^i 


M 

M c 


■gs 


-l 


where 77 = 0.1 770.1 is the radiative efficiency of accretion, K es the electron scattering (Thomson) opacity. 
We will often use accretion rate measured in units of Eddington accretion rate: 


M 

m = —-. 

fl-^Edd 


( 2 ) 


(3) 


Additional reading: There are excellent general reviews of accretion disc physics, they can be found in references 

12], [IB] BE an d [551 . 


II. DISC DRIVING MECHANISM; VISCOSITY 

In recent years there have been an impressive progress in understanding the physical mechanisms that drive disc 
accretion. It is now obvious that the turbulence in ionized Keplerian discs is due to the Magneto-Rotational Instability 
(MRI) also known as the Balbus-Hawley mechanisms [8], [7j. However, despite these developments, numerical simu¬ 
lations, even in their global, 3D form suffer still from weaknesses that make their direct application to real accretion 
flows almost infeasible. 

One of the most serious problems is the value of the ratio of the (vertically averaged) total stress to thermal 
(vertically averaged) pressure 


a = 


(P), 


(4) 


which according to most MRI simulation is ~ 10~ 3 whereas observations of dwarf nova decay from outburst un¬ 
ambiguously show that a ~ 0.1 — 0.2 51] . 2TJ|. Only recently Hirose et al. [22] showed that effects of convection 
at temperatures ~ 10 4 K increase a to values ~ 0.1. This might solve the problem of discrepancies between the 
MRI-calculated and the observed value of a HU. One has, however, to keep in mind that the simulations in question 
have been performed in a so-called shearing box and their validity in a generic 3D case has yet to be demonstrated. 

Another problem is related to cold discs such as quiescent dwarf nova discs HU or protostellar discs [T]. For the 
standard MRI to work, the degree of ionization in a weakly magnetized, quasi-Keplerian disc must be sufficiently high 
to produce the instability that leads to a breakdown of laminar flow into turbulence which is the source of viscosity 
driving accretion onto the central body. In cold discs the ionized fraction is very small and might be insufficient for 
the MRI to operate. In any case in such a disc non-ideal MHD effects are always important. All these problems still 
await their solution. 

Finally, and very relevant to the subject of this Chapter there is the question of stability of discs in which the 
pressure is due to radiation and opacity to electron scattering. According to theory such discs should be violently 
(thermally) unstable but observations of systems presumed to be in this regime totally infirm this prediction. MRI 
simulations not only do not solve this contradiction but rather reinforce it [24] . 


A. The q— prescript ion 


The a-prescription [53] is a rather simplistic description of the accretion disc physics but before one is offered better 
and physically more reliable options its simplicity makes it the best possible choice and has been the main source of 
progress in describing accretion discs in various astrophysical contexts. 

One keeps in mind that the accretion-driving viscosity is of magnetic origin, but one uses an effective hydrodynam- 
ical description of the accretion flow. The hydrodynamical stress tensor is (see e.g. [34| 1 


dv u 


dfl 


Tr *- pv OR -PdlnR 1 


( 5 ) 


where p is the density, v the kinematic viscosity coefficient and v v the azimuthal velocity (v v = Rfl). 
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In 1973 Shakura & Sunyaev proposed the (now famous) prescription 

T rl p — aP, 

where P is the total thermal pressure and a < 1. This leads to 


dfi n_1 
din R 

where c s = \[Pjp is the isothermal sound speed and p the density. For the Keplerian angular velocity 

1/2 


this becomes 




v = 3 otcr s /Q, K . 


( 6 ) 


(7) 


( 8 ) 


(9) 


Using the approximate hydrostatic equilibrium Eq. (191 one can write this as 

v ~ | ac s H. (10) 

Multiplying the rhs of Eq. § by the ring length ( 2irR ) and averaging over the (total) disc height one obtains the 
expression for the total torque 


T = 2n RY,vR 


dfl 

d In R 1 


where 


For a Keplerian disc 


/• + 00 


E = 


pdz. 


T = 3ttEu£ k , 


( 11 ) 


( 12 ) 


(13) 


(£k = R 2 Qk is the Keplerian specific angular momentum.) 

The viscous heating is proportional to to T rtp (dSl /dR) (34j . In particular the viscous heating rate per unit volume is 


q + = ~T r , 


dfl 


' rv d\uR' 

which for a Keplerian disc, using Eq. ([6]), can be written as 

q+ = ^afl K P, 

and the viscous heating rate per unit surface is therefore 

Q+ = — = 9 E^. 

(The denominator in the first rhs is 2 x 2nR taking into account the existence of two disc surfaces.) 
Additional reading: References m , 0 and [22] , 


(14) 


(15) 


(16) 


III. GEOMETRICALLY THIN KEPLERIAN DISCS 

The 2D structure of geometrically thin, non-self-gravitating, axially symmetric accretion discs can be split into a 
1+1 structure corresponding to a hydrostatic vertical configuration and radial quasi-Keplerian viscous flow. The two 
ID structures are coupled through the viscosity mechanism transporting angular momentum and providing the local 
release of gravitational energy. 
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A. Disc vertical structure 


The vertical structure can be treated as a one-dimensional star with two essential differences: 

1 . the energy sources are distributed over the whole height of the disc, while in a star there limited to the nucleus, 

2 . the gravitational acceleration increases with height because it is given by the tidal gravity of the accretor, while 
in stars it decreases as the inverse square of the distance from the center. 

Taking these differences into account the standard stellar structure equations (see e.g. [48]) adapted to the descrip¬ 
tion of the disc vertical structure are listed below. 


• Hydrostatic equilibrium 

The gravity force is counteracted by the force produced by the pressure gradient: 

dP 

~r = P9z . 
dz 

where g z is the vertical component (tidal) of the accreting body gravitational acceleration: 


d \ GM 
~dz [{R 2 +Z 2 ) 1 / 2 


GM z 

~Wr' 


(17) 


(18) 


The second equality follows from the assumption that z « R. Denoting the typical (pressure or density) scale- 
height by H the condition of geometrical thinness of the disc is H/R <C 1 and writing dP/dz ~ P/H, Eq. 
can be written as 


H _ cs_ 
R VK 


(19) 


where vk 


\JGM/R is the Keplerian velocity and we made use of Eq. (18). From Eq. (19) it follows that 


H 1 

~ 77 — : 

C s It, K 


( 20 ) 


where fdyn is the dynamical time. 


• Mass conservation 

In ID hydrostatic equilibrium the mass conservation equation takes the simple form of 


dg 

dz 


= 2 p. 


( 21 ) 


• Energy transfer - temperature gradient 


For radiative energy transport 


dhiT ^dlnP 

dz dz 


_ k-rPF z 
4P r cg z ’ 


( 22 ) 


(23) 


where P r is the radiation pressure and kr the Rosseland mean opacity. From Eqs. (221 and (231 one recovers 
the familiar expression for the radiative flux 


16 crT 3 dT 4cr dT 4 


F : 


3 krp dz 3kr p dz 


(24) 
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(F z is positive because the temperature decreases with z so dT/dz < 0). 
The photosphere is at optical thickness r ~ 2/3 (see Eq. 


T = Tr, 


74). 


The boundary conditions are: z = 0, F z = 0, 
? = 0 at the disc midplane; at the disc photosphere c = E and T 4 (t = 2/3) = T 4 ff . For a detailed 


discussion of radiative transfer, temperature stratification and boundary conditions see Sect. MB 


In the same spirit as Eq. (19) one can write Eq. (24) as 


4 crT c 4 _ 8 crT c 4 
3 krpH 3 «rE’ 


(25) 


where T c is the mid-plane (“central”) disk temperature. Using the optical depth r = n^pH = (1/2)krE, this 
can be written as 


FA H) 




(26) 


(see Eq. 77 for a rigorous derivation of this formula). 


Remark 1. In some references (for example in [ 18| ') the numerical factor on the rhs is “4/3” instead of “8/3”. 
This is due to a different definition of E: in our case it is = 2 pH, whereas in [T8j E = pH. // 


In the case of convective energy transport V = V conv . Because convection in discs is still not well understood 
(see, however, [25]) there is no obvious choice for V con v In practice a prescription designed by Paczynski [5T] 
for extended stellar envelope is used m but this most probably does not represent very accurately what is 
happening in convective accretion discs HU 


• Energy conservation 

Vertical energy conservation should have the form 


dFz 

dz 


= Q + (z), 


(27) 


where q + {z) corresponds to viscous energy dissipation per unit volume. 

Remark 2. In contrast with accretion discs, stellar envelopes have dF z /dz = 0 

The a prescription does not allow deducing the viscous dissipation stratification (z dependence), it just says 
that the vertically averaged viscous torque is proportional to pressure. Most often one assumes therefore that 


Q + (z) = o^Kf’fz), 


(28) 


by analogy with Eq. (15) but such an assumption is chosen because of its simplicity and not because of some 


physical motivation. In fact MRI numerical simulations suggest that dissipation is not stratified in the way as 
pressure [551 . 


The vertical structure equations have to be completed by the equation of state (EOS): 


P = P r + P g = ^T 4 + - P T, (29) 

3c p 

where 7 Z is the gas constant and p the mean molecular weight, and an equation describing the mean opacity 
dependence on density and temperature. 


B. Disc radial structure 


• Continuity (mass conservation) equation has the form 


9E 

~dt 


IA 

RdR 


(RTlVr) + 


S(R,t) 

2nR 


(30) 


where S(R,t) is the matter source (sink) term. 
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In the case of an accretion disc in a binary system 


S(R,t) = 


dM ext (R, t) 
dR 


( 31 ) 


represents the matter brought to the disc from the Roche lobe filling/mass losing (secondary) companion of the 
accreting object. M ex t ~ where M tr is the mass transfer rate from the companion star. Most often one 
assumes that the transfer stream delivers the matter exactly at the outer disc edge, but although this assumption 
simplifies calculations it is contradicted by observations that suggest that the stream overflows the disc surface(s). 

Angular momentum conservation 


5E7 1 d , 1 d dn 

dt ~ R8R^ R Vl ^ + RdR\ R ^ U dR 


S t {R,t) 
2t tR ' 


(32) 


This conservation equation reflects the fact that angular momentum is transported through the disc by a viscous 
stress T rlf , = RYivdVl/dR. Therefore, if the disc is not considered infinite (recommended in application to real 
processes and systems) there must be somewhere a sink of this transported angular momentum Sg(R,t). For 
binary semi-detached binary systems there is both a source (angular momentum brought in by the mass transfer 
stream form the stellar companion) and a sink (tidal interaction taking angular momentum back to the orbit). 
The two respective terms in the angular momentum equation can be written as 


Sj(R,t) = 


4 dM ext T tid (7?) 


2t tR dR 


2t tR 


(33) 


Assuming = f2^, from Eqs. (301 and (32) one can obtain an diffusion equation for the surface density E: 


— = — — I Ft 1 / 2 — IvYR 1 / 2 

dt R dR 1 R dRl^ 


(34) 


Comparing with Eqs. (30) one sees that the radial velocity induced by the viscous torque is 

3 d 


ER!/2 dR L 


/E R 1/2 


which is an example of the general relation 


v 

R' 


(35) 


(36) 


Using Eq.(10) one can write 


R 


v v 


74 

v 


-i H (H 
a c s \R 


(37) 


The relation between the viscous and the dynamical times is 

-2 


4 


H\ 

R I tdyn ' 


(38) 


In thin (77/7? <C 1) accretion discs the viscous time is much longer that the dynamical time. In other words, 
during viscous processes the vertical disc structure can be considered to be in hydrostatic equilibrium. 

Energy conservation 


The general form of energy conservation (thermal) equation can be written as: 


ds ds ds + 

pT dn ■ =f,T dt +Vr lm~ q “ 


di 


(39) 
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where s is the entropy density, q + and q~ are respectively the viscous and radiative energy density, and q is 
the density of external and/or radially transported energy densities. Using the first law of thermodynamics 
Tds = dU + PdV one can write 




P 


dv r 

dr 


(40) 


where U = 3ffT c //x(7 — 1). 
Vertically averaging, but taking T 
(3 = 1) one obtains 


T c , using Eq. (30) and the thermodynamical relations from Appendix (for 


dTc 3T C (RT C 1 d(Rv z ) _ n Q + -Q- Q 
dt Vl dR /icp R dR cpT, cpE’ 


(41) 


where Q + and Q~ are respectively the heating and cooling rates per unit surface. Q = Q au t + J with Q 0 ut 
corresponding to energy contributions by the mass-transfer stream and tidal torques; J(T, £) represent radial 
energy fluxes that are a more or less ad hoc addition to the 1+1 scheme to which they do not belong since it 
assumes that radial gradients ( d/dR ) of physical quantities can be neglected. 


The viscous heating rate per unit surface can be written as (see Eq. 161 


Q + — 


(42) 


while the cooling rate over unit surface (the radiative flux) is obviously 


In thermal equilibrium one has 


or = oiis. 

Q + = Q~. 


(43) 

(44) 


The cool ing time can be easily estimated from Eq. (441. The energy density to be radiated away is 
Eqs 228 and 232), so the energy per unit surface is ~ T,c 2 and the cooling (thermal) time is 


Ec? T,c 2 . , , 

*th — q _ ~ ~Q+ ~ a ^K ~ a ^ d y 


pc 2 s (see 


(45) 


Since a < 1, fth > fdyn and during thermal processes the disc can be assumed to be in (vertical) hydrostatic 
equilibrium. 

For geometrical thin (H/R -C 1) accretion discs one has the following hierarchy of the characteristic times 


^dvn + fth ^ t\ 


(46) 


(This hierarchy is similar to that of characteristic times in stars: the dynamical is shorter than the thermal 
(Kelvin-Helmholtz) and the thermal is much shorter than the thermonuclear time.) 


C. Self-gravity 


In this Chapter we are interested in discs that are not self-gravitating, i.e. in discs where the vertical hydrostatic 
equilibrium is maintained against the pull of the accreting body’s tidal gravity whereas the disc’s self-gravity can be 
neglected. We will see now under what conditions this assumption is satisfied. 

The equation of vertical hydrostatic equilibrium can be written as 


IdP 

-~r = -g = \~9z - g s ) = g z 

p dz 


i + ^ 


— : ~ 9z (1 + A ), 


(47) 


therefore self-gravity is negligible when A 1. Treating the disc as an infinite uniform plane (i.e. assuming the surface 
density does not vary too much with radius) one can write its self gravity as g s = 27rG£, whereas the ^-component 


of the gravity provided by the central body is g z = VL 2 k z (Eq. 18). Therefore evaluating A at z = H one gets 


A h ■■= ~ 


_ 2ttG£ 

H & 2 k H 


(48) 













Ah is related to the so-called Toomre parameter [561 


Qt 


c s n 

nGY, ’ 


(49) 


widely used in the studies of gravitational stability of rotating systems, through Ah ~ Q T . We will therefore express 


the condition of negligible self-gravity (gravitational stability) as 

Qt > 1- 


Using Eqs. (19), (10) and (57) one can write the Toomre parameter as 


Qt = 


or as function of the mid-plane temperature T = 10 4 X4 


GM 1 


Qt « 4.6 x 10 7 


aT: 


,3/2 


TO TO 


(50) 


(51) 


(52) 


where m = M/Mq. This shows that hot ionized (T > 10 4 ) discs become self-gravitating for high accretor masses 
and high accretion rates. Discs in close binary systems (to < 30) are never self-gravitating for realistic accretion 
rates (to < 1000, say) and even in IMBH binaries (if they exist) (hot) discs would also be free of gravitational 
instability. Around a supermassive black hole, however, discs can become self-gravitating quite close to the black 
hole. For example when the black hole mass is m = 10 8 a hot disc will become self-gravitating at R/Rs ~ 100, for 
to ~ 1(R 2 . In general, geometrically thin, non-self-gravitating accretion discs around supermassive black holes have 
very a limited radial extent. 

Additional reading: References !12), 19;.[20 ■ [35], [43] and [56] . 


D. Stationary discs 


In the case of stationary (d/dt = 0) discs Eq. (30) can be easily integrated giving 


M := 2nRT,v r , 

where the integration constant M (mass/time) is the accretion rate. 
Also the angular momentum equation (32) can be integrated to give 


Or, using Eq. (531 


— 27t RT,v r l + 2 ttR?YiI'Q! = const. 


— M£ + T = const., 


(53) 


(54) 


(55) 


where the torque T := 27ri? 3 E^^l , ; (for a Keplerian disc T = 37ri? 2 Ei/fljf). 

Assuming that at the inner disc radius the torque vanishes one gets const. = —Ml ; n , where l\ n is the specific 
angular momentum at the disc inner edge. Therefore 


M(l - 4.) = T 

which is a simple expression of angular momentum conservation. 

For Keplerian discs one obtains an important relation between viscosity and accretion rate 


V ~ Ak 


1 - 


Rin 

R 


1 / 2 ' 


From Eqs. (57), (42), (43), and the thermal equilibrium equation (44) it follows that 

1/2I 


= x- 


8 crTy 


3 r 


3 GM M 
87t R 3 


1 - 


Rm 

R 


(56) 


(57) 


(58) 
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FIG. 1: The observed temperature profile of the accretion disc of the dwarf nova Z Cha in outburst. Near the outburst maximum 
such a disc is in quasi-equilibrium. The observed profile, represented by dots (pixels), is compared with the theoretical profiles 
calculated from Eq. (581 and represented by continuous lines. Pixels with R < 0.03-Rli correspond to to the surface of the 
accreting white dwarf whose temperature is 40 000 K. The accretion rate in the disc is « 10 -9 M@y -1 . [Figure 6 from |23|]. 


This relation assumes only that the disc is Keplerian and in thermal ( Q + = Q~) and viscous (M = const.) 
equilibrium. The viscosity coefficient is absent because of the thermal equilibrium assumption: in such a state the 
emitted radiation flux cannot contain information about the heating mechanism, it only says that such mechanism 
exists. Steady discs do not provide information about the viscosity operating in discs or the viscosity parameter a. 
To get this information one must consider (and observe) time-dependent states of accretion discs. 

From Eq. (581 one obtains a universal temperature profile for stationary Keplerian accretion discs 


T eS ~ ir 3 / 4 . 


(59) 


For an optically thick disc the observed temperature T ~ T e g and T ~ R~ 3 / 4 should be observed if stationary, 
optically thick Keplerian discs exist in the Universe. And vice versa, if they are observed, this proves that such discs 
exist not only on paper. In 1985 Horne & and Cook ([23! presented the observational proof of existence of Keplerian 
discs when they observed the dwarf nova binary system ZCha during outburst (see Fig. [I]). 


1. Total luminosity 

The total luminosity of a stationary, geometrically thin accretion disc, i.e. the sum of luminosities of its two surfaces, 


is 


r R o 


crT 4 ff 2 t xRdR = 


3 GMM r R ° M 


'Rir 


1 - 


Rip 

It 


1 / 2 ' 


dR 

W 


(60) 


For R out —► oo this becomes 


L disc — 2 R 


1 GMM 1 


= rK 


(61) 


In the disc the radiating particles move on Keplerian orbits hence they retain half of the potential energy. If the 
accreting body is a black hole this leftover energy will be lost (in this case, however, the non-relativistic formula of 
Eq. 61 does not apply - see Eq. 178 ) In all the other cases the leftover energy will be released in the boundary layer, 
if any, and at the surface of the accretor, from where it will be radiated away. 
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Radius/R u 


FIG. 2: The observed temperature profile of the accretion disc of the dwarf nova Z Cha in quiescence. This one of the most 
misunderstood figures in astrophysics (see text). In quiescence the disc in not in equilibrium. The flat temperature profile is 
exactly what the disc instability model predicts: in quiescence the disc temperature must be everywhere lower than the critical 
temperature, but this temperature is almost independent of the radius (see Eq. 


891 . [Figure 11 from [58 


The factor “3” in the rlrs of Eq. (58) shows that radiation by a given ring in the accretion disc does not come only 

(62) 


from local energy release. Indeed, in a ring between R and R + dR only 

GMMdR 


is being released, while 


2 x 2nRQ + dR = 


2 R? 


3 GMM 
2 R 2 


1 - 


Rin 

R 


1/2 


dR 


is the total energy release. Therefore the rest 

GMM 


R 2 


^ _ 3 / R in 
2 l R 


1 / 2 ' 


dR 


(63) 


(64) 


must diffuse out from smaller radii. This shows that viscous energy transport redistributes energy release in the disc. 


E. Radiative structure 


Here we will show an example of the solution for the vertical thin disc structure which exhibit properties impossible 
to identify when the structure is vertically averaged. We will also consider here an irradiated disc - such discs are 
present in X-ray sources. 

We write the energy conservation as : 

dF , 

— =q + {R,z), (65) 


where F is the vertical (in the z direction) radiative flux and q + {R,z) is the viscous heating rate per unit volume. 
Eq. ( |65| states that an accretion disc is not in radiative equilibrium ( dF/dz ^ 0), contrary to a stellar at mosph ere. 
For this equation to be solved, the function q + (R,z) must be known. As explained and discussed in Sect, 
viscous dissipation is often written as 


III A the 


q + (R, z) = ^aCl K P(z). 


( 66 ) 
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Viscous heating of this form has important implications for the structure of optically thin layers of accretion discs 
and may lead to the creation of coronae and winds. In reality it is an an hoc formula inspired by Eq. ©• We don’t 
know yet (see, however, ED how to describe the viscous heating stratification in an accretion disc and Eq. ( |66| ) just 
assumes that it is proportional to pressure. It is simple and convenient but it is not necessarily true. 

When integrated over z, the rhs of Eq. (651 using Eq. ( 66 ) is equal to viscous dissipation per unit surface: 


o r+oo 

F + = -afl K / Pdz, 
Jo 


(67) 


where F + = (1/2 )Q + because of the integration from 0 to +oo while Q + contains E which is integrated from —oo to 
+oo (Eq. [T2|. 


One can rewrite Eq. (65) as 


dF F ■ 

~r = -/O)—: 

dr r to t 


( 68 ) 


where we introduced a new variable, the optical depth dr = —nnpdz, kr, being the Rosseland mean opacity and 
r to t = f 0 +X Knpdz is the total optical depth. f(r) is given by: 


/ 0 ) = 


P 


(fo + °° K RPdz) 


(sr?*-- 


KR P 


(69) 


As p decreases approximately exponentially, /(r) is the ratio of two rather well defined scale heights, the pressure 
and the opacity scale heights, which are comparable, so that / is of order of unity. 

At the disc midplane, by symmetry, the flux must vanish: J 7 '(rtot) = 0, whereas at the surface, (r = 0) 


F(0) = oT* s = F+. 


(70) 


Equation |70[ ) states that the total flux at the surface is equal to the energy dissipated by viscosity (per unit time 
and unit surface). The solution of Eq. ( 68 ) is thus 


F(t) = F+ 1 - 


fo f( T )dT 

7"tot 


(71) 


where JJ tot f(r)dT = r to t- Given that / is of order of unity, putting /(r) = 1 is a reasonable approximation. The 
precise form of f(r) is more complex, and is given by the functional dependence of the opacities on density and 
temperature; it is of no importance in this example. We thus take: 


F{t) =F+ ( 1- 

Etot 


(72) 


To obtain the temperature stratification one has to solve the transfer equation. Here we use the diffusion approxi¬ 
mation 




(73) 


appropriate for the optically thick discs we are dealing with. The integration of Eq. (73) is straightforward and gives 


T\r) - T‘( 0) = | T D - 2^ I 


(74) 


The upper (surface) boundary condition is: 


T 4 (0) = -T e 4 ff + 7T 


4 

irr ’ 


(75) 


where 7] 4 r is the irradiation temperature, which depends on r, the albedo, the height at which the energy is deposited 
and on the shape of the disc. In Eq. (75) T(0) corresponds to the emergent flux and, as mentioned above, T e g 
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corresponds to the total flux (crT 4 ff = 
is thus : 


Q + ) which explains the factor 1/2 in Eq (75). The temperature stratification 




2 

3 


+ T- 

I ir 


(76) 


For r to t 1 the first term on the rhs has the form familiar from the stellar atmosphere models in the Eddington 
approximation. 

In this case at r = 2/3 one has T{ 2/3) = T e g. 

Also for rtot /§> 1, the temperature at the disc midplane is 

T c 4 = T 4 (r tot ) = ^r tot T 4 ff + I) 4 r - (77) 

It is clear, therefore, that for the disc inner structure to be dominated by irradiation and the disc to be isothermal 
one must have 

p. 

— = —^ » F+ (78) 

Ttot 7”tot 

and not just El rr F + as is sometimes assumed. The difference between the two criteria is important in LMXBs 
since, for parameters of interest, rtot > 10 2 — 10 3 in the outer disc regions. 


F. Shakura-Sunyaev solution 


In their seminal and famous paper Shakura & Sunyaev |53j , found power-law stationary solutions of the simplified 
version of the thin-disc equations presented in Sects. |III A IIIB and HID The 8 equations for the 8 unknowns T c , 
p , P, E, H , v, t and c s can be written as 


S = 2 Hp 
H= CsR3/2 


(GM) 1 /2 


(i) 

(n) 



(m) 


P = + l^T 4 

p, 3c 

t(p,S,T c ) = K R (p,T c ) E 


i/(p,E,T c ,a) = 7 otc s H 


(iv) 

(v) 

(vi) 


v M 
" S= 3^ 




1 / 2 ' 


(vn) 


8 oTl 
3 T 


3 GMM 
8tt R 3 


‘-'I 


1 / 2 ' 


(vm) 


Equations (i) and (n) correspond to vertical structure equations (21) and (19), Eq. (vn) is the radial Eq. (57), while 
Eq. (vm) connects vertical to radial equations. Eq. (m) defines the sound speed, Eq. (iv) is the equation of state and 
(vi) contains the information about opacities. The viscosity a parametrization introduced in [53j provides the closure 
of the 8 disc equations. Therefore they can be solved for a given set of a, M, R and M. 

Power-law solutions of these equations exist in physical regimes where the opacity can be represented in the Kramers 
form k = K 0 p n T m and one of the two pressures, gas or radiation, dominates over the other. In [53] three regimes have 
been considered: 
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a. ) P T P g and K es » Kg 

b. ) P g ^> P r and K es Kg 

c. ) P g P r and Kg K es . 

Regimes a.) and b.) in which opacity is dominated by electron scattering will be discussed in Sect. [V] Here we will 


present the solutions of regime c.), i.e. we will assume that 

P r = 0 and kr = ng = 5 x 10 2 A pT ~ 7 ^ 2 cm 2 g -1 . (79) 

The solution for the surface density E, central temperature T c and the disc relative height (aspect ratio) are respectively 

E = 23a- 4/5 m 1 / 4 i?” 3/4 M 1 7 7 /10 / 7/10 gcm- 2 , (80) 

T c = 5.8 x 10 4 a~ 1/5 m 1/A R^ /4 Mi( 10 f 3/10 K, (81) 

^ = 2.4 x 10“ 2 a _ 1 / 10 m _ 3 / 8 i?^ 8 M 1 3 7 /20 / 3/2 °, (82) 


where m = M/Mq, i?io = R/(10 10 cm), M 17 = M/(10 17 gs 4 ), and / = 1 — (Rin/R) 1 / 2 . 



FIG. 3: Stationary accretion disc surface density profiles for 4 values of accretion rate. From top to bottom: M = 10 18 ,10,10 16 
and 10 15 gs _1 . m = IOMq, a = 0.1. The continuous line corresponds to the un-irradiated disc, the dotted lines to an irradiated 


configuration. The inner, decreasing segments of the continuous lines correspond to Eq. (80). Dashed lines describe irradiated 
disc equilibria (see Sect. IV C|) [Figure 9 from llfil l. 


Although for a 10 M 0 black hole, say, Shakura-Sunyaev solutions (80), (81) and (81) describe discs rather far from 
its surface (R > 10 4 Rs) the regime of physical parameters it addresses, especially temperatures around 10 4 K are of 
great importance for the disc physics because it is where accretion discs become thermally and viscously unstable. 
This instability triggers dwarf nova outbursts when the accreting compact object is a white dwarf and (soft) X-ray 
transients in the case of accreting neutron stars and black holes. 

It is characteristic of the Shakura-Sunyaev solution in this regime that the three E, T c and T e g radial profiles vary 
as i? -3 / 4 . (This implies that the optical depth r is constant with radius - see Eq. vm.) For high accretion rates 
and small radii the assumption of opacity dominated by free-free and bound-free absorption will break down and the 
solution will cease to be valid. We will come to that later. Now we will consider the other disc end: large radii. 

One sees in Fig. [5J that for given stationary solution (M = const.) the A -3 / 4 slope of the E profiles extends down 
only to a minimum value E m i n (A) after which the surface density starts to increase. With the temperature dropping 
below 10 4 K the disc plasma recombines and there is a drastic change in opacities leading to a thermal instability. 

Additional reading: We have assumed that accretion discs are flat. This might not be true in general because 
accretion discs might be warped. This has important and sometimes unexpected consequences; see e.g [ 28] [40] and 
1451. and references therein. 
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IV. DISC INSTABILITIES 

In this section we will present and discuss the disc thermal and the (related) viscous instabilities. First we will 
discuss in some detail the cause of the thermal instability due to recombination. 


A. The thermal instability 


A disc is thermally stable if radiative cooling varies faster with temperature than viscous heating. In other words 


d In crT^ff d In Q + 
d In T r ^ d In T r 


(83) 


Using Eq. (771 one obtains 


dlnT c % 

dlnT c 


= 4 


1 - 


T irr 

T r 


-l 


din k 
dlnT c 


(84) 


In a gas pressure dominated disc Q + ~ pT H ~ E T ~ T c . The thermal instability is due to a rapid change of 
opacities with temperature when hydrogen begins to recombine. At high temperatures d In n/d In T c ss —4 (see Eq. 
79). In the instability region, the temperature exponent becomes large and positive din At/d In T c ss 7— 10, and in the 
end cooling is decreasing with temperature. One can also see that irradiation by furnishing additional heat to the 
disc can stabilize an otherwise unstable equilibrium solution (dashed lines in Fig. [3|. 

This thermal instability is at the origin of outbursts observed in discs around black-holes, neutron stars and white 
dwarfs. Systems containing the first two classes of objects are known as Soft X-ray transients (SXTs, where “soft” 
relates to their X-ray spectrum), while those containing white-dwarfs are called dwarf-novae (despite the name that 
could suggest otherwise, nova and supernova outbursts have nothing to do with accretion disc outbursts). 


B. Thermal equilibria: the S-curve 


We will first consider thermal equilibria of an accretio n disc in which heating is due only to local turbulence, 
leaving the discussion of the effects of irradiation to Section 


IV C 


We put therefore Ti rr = Q = 0. Such an assumption 
corresponds to discs in cataclysmic variables which are the best testbed for standard accretion disc models. The 


thermal equilibrium in the disc is defined by the equation Q = Q + (see Eq. 41), i.e. by 


crT cS — -z/EO K 


(85) 



2 (g cm 2 ) 

FIG. 4: Thermal equilibria of a ring in an accretion discs around a m = 1.2 white dwarf. The distance from the center is 
10 9 cm; accretion rate 6.66 x 10 16 g/s. The solid line corresponds to Q + = Q~. £ m i n is the critical (minimum) surface density 
for a hot stable equilibrium; E max the maximum surface density of a stable cold equilibrium. 
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(Eq. |16[ ). In general, v is a function of density and temperature and in the following we will use the standard 
a-prescription Eq. ([7| . The energy transfer equation provides a relation between the effective and the disc midplane 
temperatures so that thermal equilibria can be represented as a T e g (E) - relation (or equivalently a M(E)-relation). 
In the temperature range of interest (10 3 < T e g < 10 5 ) this relation forms an 5 on the (E, T e g ) plane as in Fig. [4j The 
upper, hot branch corresponds to the Shakura-Sunyaev solution presented in Section |III F| The two other branches 
correspond to solutions for cold discs - along the middle branch convection plays a crucial role in the energy transfer. 

Each point on the (E,T e fj) 5-curve represents an accretion disc’s thermal equilibrium at a given radius, i.e. a 
thermal equilibrium of a ring at radius R. In other words each point of the 5-curve is a solution of the Q + = Q~ 
equation. Points not on the 5-curve correspond to solutions of Eq. (41) out of thermal equilibrium: on the left of 
the equilibrium curve cooling dominates over heating, Q + < Q~\ on the right heating over cooling Q + > Q~. It 
is easy to see that a positive slope of the T e ff(E) curve corresponds to stable solutions. Indeed, a small increase of 
temperature of an equilibrium state (an upward perturbation) on the upper branch, say, will bring the ring to a state 
where Q + < Q~ so it will cool down getting back to equilibrium. In a similar way an downward perturbation will 
provoke increased heating bringing back the system to equilibrium. 

The opposite is happening along the 5-curve’s segment with negative slope as both temperature increase and 
decrease lead to a runaway. The middle branch of the 5-curve corresponds therefore to thermally unstable equilibria. 

A stable disc equilibrium can be represented only by a point on the lower, cold or the upper, hot branch of the 
5-curve. This means that the surface density in a stable cold state must be lower than the maximal value on the 
cold branch: E max , whereas the surface density in the hot stable state must be larger than the minimum value on 
this branch: E m j n . Both these critical densities are functions of the viscosity parameter a, the mass of the accreting 
object, the distance from the center and depend on the disc’s chemical composition. In the case of solar composition 
the critical surface densities are 


S min (5) = 39.9 c^' 80 Rio 11 m~ 0 ' 37 gcm -2 (86) 

£ max (5) = 74.6 c^- 83 R\ 0 18 mi 0A0 gem" 2 , (87) 

(a = O.lao.i) and the corresponding effective temperatures (T + designates the temperature at E m ; n , T~ at E max ) 

T+ = 6890 Rl Q 0 - 09 M 903 K (88) 

T~ s = 5210 i?C 0 010 Mf 04 K. (89) 


The critical effective temperatures are practically independent of the mass and radius because they characterize the 
microscopic state of disc’s matter (e.g. its ionization). On the other hand the critical accretion rates depend very 
strongly on radius: 


M+ it (i?) = 8.07 x 10 15 o^- 01 R\ o 64 Mf°' 89 gs -1 (90) 

M~ it (R) = 2.64 x 10 15 a 9 ;? 1 R w 8 Aff 9 - 85 gs" 1 . (91) 


A stationary accretion disc in which there is a ring with effective temperature contained between the critical values 
of Eq. (89) and (88) cannot be stable. Since the effective temperature and the surface density both decrease with 
radius, the stability of a disc depend on the accretion rate and the disc size (see Fig. [3^. For a given accretion rate a 
stable disc cannot have an outer radius larger than the value corresponding to Eq. ( |86[ ). 

A disc is stable if the rate (mass-transfer rate in a binary system) at which mass is brought to its outer edge 


(R ~ Rd) is larger than the critical accretion rate at this radius M A it (Rd). 

In general, the accretion rate and the disc size are determined by mechanisms and conditions that are exterior to 
the accretion process itself. In binary systems, for instance, the size of the disc is determined by the masses of the 
system’s components and its orbital period while the accretion rate in the disc is fixed by the rate at which the stellar 
companion of the accreting object loses mass, which in turn depends on the binary parameters and the evolutionary 
state of this stellar mass donor. Therefore the knowledge of the orbital period and the mass-transfer rate should 
suffice to determine if the accretion disc in a given interacting binary system is stable. Such knowledge allows testing 
the validity of the model as we will show in the next section. 


1. Dwarf nova and X-ray transient outbursts 


• Local view: the limit cycle 
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a^-O.l. <*^,,,=0.02, R=10* cm 



E (g cm-*) 


FIG. 5: Local limit cycle of the state of disc ring at 10 9 cm during a dwarf nova outbursts. The arrows show the direction of 
motion of the system in the T e ff(£) plane. The figure represents results of the disc instability model numerical simulations. 
As required by the comparison of the model with observations the values of the viscosity parameter a on the hot and cold 
branches are different. [Figure adapted from [36] ] 


Let us first describe what is happening during outbursts with a disc’s ring. Its states are represented by a point 
moving in the E — T e g plane as shown on Fig. which represents accretion disc states at R = 10 9 cm (the accreting 
body has a mass of 1.2M 0 ). To follow the states of a ring during the outburst let us start with an unstable equilibrium 
state on the middle, unstable branch and let us perturb it by increasing its temperature, i.e. let us shift it upwards 
in the T e g(E) plane. As we have already learned, points out of the S-curve correspond to solutions out of thermal 
equilibrium and in the region to the right of the S'-curve heating dominates over cooling. The resulting runaway 
temperature increase is represented by the point moving up and reaching (in a thermal time) a quasi-equilibrium 
state on the hot and stable branch. It is only a quasi-equilibrium because the equilibrium state has been assumed to 
lie on the middle branch which corresponds to a lower temperature (and lower accretion rate - see Eq. 58). Trying 
to get to its proper equilibrium the ring will cool down and move towards lower temperatures and surface densities 
along the upper equilibrium branch (in a viscous time). But the hot branch ends at E m ; n , i.e. at a temperature higher 
(and surface density lower) than required so the ring will never reach its equilibrium state. Which is not surprising 
since this state is unstable. Once more the ring will find itself out of thermal equilibrium but this time in the region 
where cooling dominates over heating. Rapid (thermal-time) cooling will bring it to the lower cool branch. There, 
the temperature is lower than required so the point representing the ring will move up towards E max where it will 
have to interrupt its (viscous-time) journey having reached the end of equilibrium states before getting to the right 
temperature. It will find itself out of equilibrium where heating dominated over cooling so it will move back to the 
upper branch. 

Locally, the state of a ring performing a limit cycle on the E -T e g plane, moves in viscous time on the stable S'-curve 
branches and in a thermal time between them when the ring is out of thermal equilibrium. The states on the hot 
branch correspond to outburst maximum and the subsequent decay whereas the quiescence correspond to moving on 
the cold branch. Since the viscosity is much larger on the hot than on the cold branch, the quiescent is much longer 
than the outburst phase. The full outburst behaviour can be understood only by following the whole disc evolution 
(Figs. [7]&||). 


C. Irradiation and black—hole X-ray transients 

We will present the global view of thermal-viscous disc outbursts for the case of X-ray transients. The main 
difference between accretion discs in dwarf novae and these systems is the X-ray irradiation of the outer disc in the 
latter. Assuming that the irradiating X-rays are emitted by a point source at the center of the system, one can write 
the irradiating flux as 

(tT\ f r = with L x = ??min (M in ,M E dd) c 2 , (92) 

where C = 1CP 3 C3, q is the radiative efficiency (which can be <C 0.1 for ADAFs - see below), M ln the accretion rate at 
the inner disc’s edge. Since the physics and geometry of X-ray self-irradiation in accreting black-black hole systems 
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FIG. 6: Example S-curves for a pure helium disk with varying irradiation temperature Tirr. The various sets of S-curves 
correspond to radii R = 10 6 , 10 u and 10 10 cm. For each radius, the irradiation temperature Ti rr is OK, 10 000K and 20 000 K. 
a = 0.16. The instable branch disappears for high irradiation temperatures. [From 211 ■ Reproduced with permission from 
Astronomy & Astrophysics, ©ESO] 


is still unknown, the best we can do is to parametrize our ignorance by q constant C that observations suggest is 
~ 10“ 3 . Of course one should keep in mind that in reality C might not be a constant m- 

Because the viscous heating is ~ M/R 3 there always exists a radius R in - for which oT > Q + = oTf. If R irr < R dl 
where R d is the outer disc radius, the o uter d isc emission will be dominated by reprocessed X-ray irradiation and 
the structure modified as shown in Sect. |IIIE Irradiation will also stabilize outer disc regions (Eq. [84] and Fig. [6]) 
allowing larger discs for a given accretion rate (see Fig. [3]). 

Irradiation modifies the critical values of the hot disc parameters: 


y+ . . 

^irr — 

rr.irr,+ _ 

1 cff ~ 


72.4 cz° 3 28 aff 8 R ii 2 Mf 0 ' 19 gcm“ 2 
2860 CZ 3 0 09 a#? 1 K[- 1 °' 15 M? m K 


M+ = 2.3 x 10 


3 u o.: 

17 y>-0.36 
«4-3 


*11 

0.04 

* 0.1 


Ri f M- 




(93) 

(94) 

(95) 


As we will see in a moment, irradiation also strongly influences the shape of outburst’s light-curve. 


* Rise to outburst maximum 

During quiescence the disc’s surface density, temperature and accretion rate are everywhere (at all radii) on the cold 
branch, below their respective critical values E max (i?), T~ s and M~ it (R). It is important to realize that in quiescence 
the disc is not steady: Ad ^ const. Matter transferred from the stellar companion accumulates in the disc and is 
redistributed by viscosity. The surface density and temperature increase (locally, this means that the solution moves 
up along the lower branch of the S-curve) finally reaching their critical values. In Fig. [ 7 ] this happens at ~ 10 10 cm. 
The disc parameters entering the unstable regime triggers an outburst. In the local picture this corresponds to leaving 
the lower branch of the S'-curve. The next ‘moment’ (in a thermal time is represented in the left panels of Figure [ 7 ] 
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IV C 


The upper left panel shows M- m and Mi rr (dotted line); the bottom left 
Each dot corresponds to one of the E and T c profiles in the right panels. The heating front 


FIG. 7: The rise to outburst described in Sect, 
panel shows the V magnitude. 

propagates outwards. The disc expands during the outburst due to the angular momentum transport of the material being 
accreted. At t ~ 5.5 days the thin disc reaches the minimum inner disc radius of the model. The profiles close to the peak 
are those of a steady-state disc (E oc T c oc R~ 3 ' 4 ). [From (T5]. Reproduced with permission from Astronomy & Astrophysics, 
©ESO] 


This is when a large contrast forms in the midplane temperature profile and when a surface-density spike is already 
above the critical line. The disc is undergoing a thermal runaway at r ~ 8 x 10 9 cm. The midplane temperature rises 
to ~ 70000 K. This raises the viscosity which leads to an increase of the surface-density and a heating fronts start 
propagating inwards and outwards in the disc as seen in Fig. [7J In this model the disc is truncated at an inner radius 
i?i n w6x 10 9 cm so the inwards propagating front quickly reaches the inner disc radius with no observable effects. It 
is the outwards propagating heating front that produces the outburst by heating up the disc and redistributing the 
mass and increasing the surface density behind it because it is also a compression front. 

One should stress here that two ad hoc elements must be added to the model for it to reproduce observed outbursts 
of dwarf novae and X-ray transients. 

• Viscosity. First, if the increase in viscosity were due only to the rise in the temperature through the speed 
of sound (y oc c 2 , see Eq. |9| the resulting outbursts would have nothing to do with the observed ones. To 
reproduce observed outbursts one increases the value of a when a given ring of the disc gets to the hot branch. 
Ratios of hot-to-cold a of the order of 4 are used to describe dwarf nova outburst. Although in the outburst 
model the a increase is an ad hoc assumption, recent MRI simulations with physical parameters corresponding 
to dwarf nova discs show an a increase induced by the appearance of convection |22l . 

• Inner truncation. Second, as mentioned already, the inner disc is assumed to be truncated in quiescence and 
during the rise to outburst. Although such truncation is implied and/or required by observations, its physical 
origin is still uncertain. The inner part of the accretion flow is of course not empty but supposed to form a 
M = const. ADAF (see Sect. [v| . 

In our case (Fig. [?|, the heating front reaches the outer disc radius. This corresponds to the largest outbursts. 
Smaller-amplitude outbursts are produced when the front does not reach the outer disc regions. In an inside-out 
outburst 1 the surface-density spike has to propagate uphill, against the surface-density gradient because just before 
the outburst E ~ R 118 - roughly parallel to the critical surface-density. Most of the mass is therefore contained in 
the outer disc regions. A heating front will be able to propagate if the post-front surface-density is larger than X m ; n 
in other words, if it can bring successive rings of matter to the upper branch of the S'-curve. If not, a cooling front 
will appear just behind the E spike, the heating front will die-out and the cooling front will start to propagate inwards 
(the heating-front will be ‘reflected’). 


1 X-ray transient outbursts are always of inside-out type. In dwarf novae both inside-out and outside-in outbursts are observed and result 
from calculations [311 . 
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FIG. 8: Decay from outburst peak. The decay is controlled by irradiation until evaporation sets in at t ~ 170 days (Mi n = 
Mevap(-Rmin)). This cuts off irradiation and the disc cools quickly. The irradiation cutoff happens before the cooling front 
can propagate through most of the disc, hence the irradiation-controlled linear decay (t ~ 80 — 170 days) is not very visible 
in the lightcurve. Ti rr (dotted line) is shown for the last temperature profile. [From |15| . Reproduced with permission from 
Astronomy & Astrophysics, ©ESO] 


The difficulty inside-out fronts encounter when propagating is due to angular-momentum conservation. In order 
to move outwards the X-spike has to take with it some angular momentum because the disc’s angular momentum 
increases with radius. For this reason inside-out front propagation induces a strong outflow. In order for matter to 
be accreted, a lot of it must be sent outwards. That is why during an inside-out dwarf-nova outburst only ~ 10% of 
the disc’s mass is accreted onto the white dwarf. In X-ray transients irradiation facilitates heating front propagation 
(and disc emptying during decay - see next section). 

The arrival of the heating front at the outer disc rim does not end the rise to maximum. After the whole disc is 
brought to the hot state, a surface density (and accretion rate) ‘excess’ forms in the outer disc. The accretion rate 
in the inner disc corresponds to the critical one but is much higher near the outer edge. While irradiation keeps the 
disc hot the excess diffuses inwards until the accretion rate is roughly constant. During this last phase of the rise to 
outburst maximum Afj n increases by a factor of 3: 

Mm ax « 3Mj+ « 7.0 x 10 17 CI 3 0 - 36 R 2 d jl m " 064 g s " 1 . (96) 

Irradiation has little influence on the actual vertical structure in this region and T c oc X oc as in a non- 

irradiated steady disc. Only in the outermost disc regions does the vertical structure becomes irradiation-dominated, 
i.e. isothermal. 


* Decay 

Fig. H shows the sequel to what was described in Fig. [7j In general the decay from the outburst peak of an 
irradiated disc can be divided into three parts: 

• First, X-ray irradiation of the outer disc inhibits cooling-front propagation. But since the peak accretion rate 
is much higher than the mass-transfer rate , 2 the disc is drained by viscous accretion of matter. 

• Second, the accretion rate becomes too low for the X-ray irradiation to prevent the cooling front from propa¬ 
gating. The propagation speed of this front, however, is controlled by irradiation. 

• Third, irradiation plays no role and the cooling front switches off the outburst on a local thermal time-scale. 


2 The peak luminosity is ~ 3A7© ( R d ): and the for the disc to be unstable the mass-transfer rate must be lower than the critical rate: 

Mu < M+(R d ). 
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‘Exponential decay’ 


In Fig. [8] the “exponential decay” the phase lasts until roughly day 80-100. At the outburst peak the accretion rate 
is almost exactly constant with radius; the disc is quasi-stationary. The subsequent evolution is self-similar: the disc’s 
radial structure evolves through a sequence of quasi-stationary (M(r) = const) states. Therefore uT, ~ M- 1Tl (t) /3-k 
and the total mass of the disk is thus 


M d 


J 2nRY,dR oc M- m 


2 r 

3 v 


dr. 


(97) 


At the outburst peak the whole disc is wholly ionized and except for the outermost regions its structure is very 
well represented by a Shakura-Sunyaev solution. In such discs, as well as in irradiation domina ted discs, the viscosity 
coefficient satisfies the relation v oc T oc Md/( 1+ P\ In hot Shakura-Sunyaev discs (3 = 3/7 (Eq. 81), and in irradiation 
dominated discs (3 = 1/3 (Eq. 92). During the first decay phase the outer disc radius is almost constant so that using 
Eq. (97) the disc-mass evolution can be written as: 


d M d 
dt 


= —Min « M$ +f> 


(98) 


showing that M- ln evolves almost exponentially, as long as M? n can be considered as constant (i.e. over about a decade 
in M\ n . ‘Exponential’ decays in the DIM are only approximately exponential. 

The quasi-exponential decay is due to two effects: 

1. X-ray irradiation keeps the disc ionized, preventing cooling-front propagation, 

2 . tidal torques keep the outer disc radius roughly constant. 


‘Linear’ decay 


The second phase of the decay begins when a disc ring cannot remain in thermal equilibrium. Locally this corre¬ 
sponds to a fall onto the cool branch of the S'-curve. In an irradiated disc this happens when the central object does 
not produce enough X-ray flux to keep the T m -(R out ) above ~ 10 4 K. A cooling front appears and propagates down 
the disc at a speed of Uf ron t ~ o^\ 1 c s . 

In an irradiated disc, however, the transition between the hot and cold regions is set by T m because a cold branch 
exists only for T; rr < 10 4 K. In an irradiated disc a cooling front can propagate inwards only down to the radius at 

which T; rr ss 10 4 K, i.e. as far as there is a cold branch to fall onto. Thus the decay is still irradiation-controlled. 

1/2 


The hot region remains close to steady-state but its size shrinks Rhot 

Tjrr (Aftot) = COnSt). 


M ir ( (as can be seen in Eq. 


92 


with 


Thermal decay 

In the model shown in Fig. [8] irradiation is unimportant after t > 170 — 190 days because r] becomes very small for 
M; n < 10 16 g-s _1 when an ADAF forms. The cooling front thereafter propagates freely inwards, on a thermal time 
scale. In this particular case the decrease of irradiation is caused by the onset of evaporation at the inner edge which 
lowers the efficiency. In general there is always a moment at which T m becomes less than 10 4 K; evaporation just 
shortens the ‘linear’ decay phase. 


D. Maximum accretion rate and decay timescale 


Now we will see that there are two observable properties of X-ray transients that, when related one to to the other, 
provide informations and c ons traints on the physical properties of the outbursting system. The first is the maximum 


accretion rate M max (Eq. 96). The second is the decay time of the X-ray flux: as we have seen, disk irradiation 


by the central X rays traps the disk in the hot, high state, and only allows a decay of M on the hot-state viscous 
timescale. This is 


R 

3v 


( 99 ) 
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which using Eq. (|9|) gives 


Taking the critical midplane temperature T+ 


_ ( GMR ) 1 / 2 
3ac 2 

16000 K one gets for the decay timescale 


( 100 ) 


32 m 1 / 2 i?y i 2 1 a 0 J, days, 


( 101 ) 


where ao.2 = a/0.2. Eliminating R between (96) and (101) gives the accretion rate through the disk at the start of 
the outburst as 


M = 5.4 x 10 17 to - 303 (i 30 a 0 .2) 4 - 78 gs" 1 , (102) 

with t = 30^30 d. Assuming an efficiency of of 10%, the corresponding luminosity is 

L = 5.0 x 10 37 ryo.irn -3 03 (f3oao.2) 4 ’ 78 ergs -1 . (103) 


E. Comparison with observations 


1. Sub-Eddington outbursts 


The peak luminosities of most of the soft X-ray transients are sub-Eddington. Eq. (102) can be written using the 
Eddington ratio m := M / M Edd as 


to = 0A27j o .i(a o .2t 3 o) 4 - 78 m 4 ’ 03 . 


(104) 


This equation shows that the outburst peak will be sub-Eddington only if the outburst decay time is relatively short 
or the accretor (black hole) mass is high, i.e. the observed decay timescale is 

t < 50 ■qo.i 21 a^m 0 ' 84 d, (105) 


in good agreement with the compilation of X-ray transients outburst durations found in (59j . This shows that the 
standard value of efficiency 770.1 — 1, and the value ao.2 — 1 deduced from observations of dwarf novae, give the 
correct order of magnitude for the decay timescale of X-ray transients (from ss 3 days to ~ 300 days). This equation 
also implies that black hole transients should have longer decay timescales than neutron star transients, all else being 
equal. Yan and Yu [55| find that outbursts last on average ~ 2.5 x longer in black hole transients than in neutron 
star transients thus confirming this conclusion. 

For sub-Eddington outbursts Eq. (1031 gives a useful relationship between distance D, bolometric flux F and 
outburst decay time t, 


Dm P c — 1.0 m 


-1.5 




(106) 


where D = D m P c Mpc and F = 10 12 Fi 2 erg s 1 cm 2 ; F = L/AnD 2 and t = 50 1 50 d. 

Eq. (106) shows that distant (D > IMpc) X-ray sources exhibiting variability typical of soft X-ray transients 


cannot contain black holes with masses superior to stellar masses 


2. Observational tests 


Finally, one can test observationally if soft X-ray transients satisfy the necessary condition for instability M tr < 
Merit (Rd), where Merit is the critical accretion rate for either non-irradiated or irradiated discs. In Fig. [9] the critical 
accretion rates (90) and (95) for respectively non-irradiated and irradiated disc around black holes are plotted as 
M(P OI b) relation. This relation was obtained from disc-radius - orbital-separation relation Rd(a) [J2], where (from 
Kepler’s law) the orbital separation a = 3.53 x 10 1o (toi + m 2 ) 4 / 3 P^ 3 cm, where m t are the masses of the components 










22 



GRO J0422+32 
A0620-00 
GRS 1009-45 
XTE Jill8+480 
GS 1124-683 
GS 1354-64 
4U 1543-57 
XTE J1550-564 
XTE J1650-500 
GRO J1655-40 
MAXI J1659-352 
GX 339-4 
4U 1705-250 
Swift J 1753.5-0127 
GRS 1915+105 
GS 2000+25 
V 404 Cyg 


O IE 1740.7-2942 
^ GRS 1758-258 


4U 1957+115 
Cyg X-1 
LMC X-1 
LMC X-3 


Orbital period (hr) 


FIG. 9: Mass transfer rate as a function of the orbital period for SXTs with black holes. The transient and persistent sources 
have been marked with respectively filled and open symbols. The shaded grey areas indicated ‘DIM hr’ and ‘DIM non irr’ 
represent the separation between persistent (above) and transient systems (below) according to the disc instability model when, 
respectively, irradiation is taken into account and when it is neglected. The horizontal dashed line indicates the Eddington 
accretion rate for a 10M© black hole. All the upper limits on the mass transfer rate are due to lower limits on the recurrence 
time. The upper limits on the mass transfer rate of 4U 1957+115 and GS 1354-64 result from lower limits on the distance 
to the sources. The three left closed arrows do not indicate actual upper limits on the orbital period of Cyg X-1, LMC X-1 
and LMC X-3. They emphasize that the radius of any accretion disk in these three high-mass XRBs is likely to be smaller 
than the one derived from the orbital period since they likely transfer mass by a (possibly focused) stellar wind instead of fully 
developed Roche lobe overflow. In the legend, the solid horizontal line separates transient and persistent systems. (The dashed 
horizontal line stresses that the persistent nature of IE 1740.7-2942 and GRS 1758-258 is unclear.) [From [13].] 


in solar units, and Phr the orbital period in hours. Against these two critical lines the actual positions of the observed 
sources are marked. The mass transfer rate being difficult to measure, a proxy in the form of the accumulation rate 

A E 

M accum = --- (107) 

trec'O c 

has been used. A E is the energy corresponding to the integrated X-ray luminosity from during an outburst and 
t rec the recurrence time of the outbursts. One can see that all low-mass-X-ray-binary (LMXB) transients are in the 
unstable part of the figure, as they should be if the model is correct. One can also see that all black hole LMXBs 
are transient. This is not true of neutron star LMXBs. Cyg X-1 in which the stellar companion of the black hole is 
a massive star is observed to be stable but according to Fig. [9] should be transient. This is not a problem because 
in such a system matter from the high-mass companion is not transferred by Roche-lobe overflow as in LMXBs, but 
lost through a stellar wind. In this case the Rd(a) relation used in the plot is not valid - the discs in such systems 
are smaller which is marked by a left-directed arrow at the symbol marking the position of this and two other similar 
objects (LMC X-1 and LMC X-3). 

Additional reading: References na, m , m , eh, and [52] . 


V. BLACK HOLES AND ADVECTION OF ENERGY 


Until now, we have neglected advection terms in the energy and momentum equations for stationary accretion 
flows. There two regimes of parameters where this assumption is not valid, in both cases for the same reason: low 
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radiative efficiency when the time for radial motion towards the black hole is shorter than the radiative cooling time. 
Low density (low accretion rate), hot, optically thin accretion flows are poor coolers and they are one of the two 
configurations were advection instead of radiation is the dominant evacuation-of-energy (“cooling”) mechanism. Such 
optically thin flows are called ADAFs, for Advection Dominated Accretion Flows. Also advection dominated are 
high-luminosity flows accreting at high rates but they are called “slim discs” to account for their property of not 
being thin but still being described as if this were not of much importance. 

We shall start with optically thin flows. 

• ADAFs 

Advection Dominated Accretion Flows’ (ADAFs) is a term describing accretion of matter with angular momen¬ 
tum, in which radiation efficiency is very low. In their applications, ADAFs are supposed to describe inflows 
onto compact bodies, such as black holes or neutron stars; but very hot, optically thin flows are bad radiators 
in general so that, in principle, ADAFs are possible in other contexts. Of course in the vicinity of black holes 
or neutron stars, the virial (gravitational) temperature is T v ; r « 5 x 10 12 (Rs/R) K, so that in optically thin 
plasmas, at such temperatures, both the coupling between ions and electrons and the efficiency of radiation 
processes are rather feeble. In such a situation, the thermal energy released in the flow by the viscosity, which 
drives accretion by removing angular momentum, is not going to be radiated away, but will be adverted towards 
the compact body. If this compact body is a black hole, the heat will be lost forever, so that advection, in this 
case, acts as sort of a ‘global’ cooling mechanism. In the case of infall onto a neutron star, the accreting matter 
lands on the star’s surface and the (reprocessed) advected energy will be radiated away. There, advection may 
act only as a ‘local’ cooling mechanism. (One should keep in mind that, in general, advection may also be 
responsible for heating, depending on the sign of the temperature gradient - in some conditions, near the black 
hole, advection heats up electrons in a two-temperature ADAF). 

In general the role of advection in an accretion flow depends on the radiation efficiency which in turns depends 
on the microscopic state of matter and on the absence or presence of a magnetic field. If, for a given accretion 
rate, radiative cooling is not efficient, advection is necessarily dominant, assuming that a stationary solution is 
possible. 


Slim discs 

At high accretion rates, discs around black holes become dominated by radiation pressure in their inner regions, 
close to the black hole. At the same time the opacity is dominated by electron scattering. In such discs 
H/R is no longer -C 1. But this means that terms involving the radial velocity are no longer negligible since 
v r ~ ac s (H/R). In particular, the advective term in the energy conservation equation v r dS/dR (see Eq. 39) 
becomes important and finally, at super-Eddington rates, dominant. When Q + = Q adv the accretion flow is 
advection dominated and called a slim disc. 


A. Advection dominated—accretion-flow toy models 


One can illustrate the fundamental properties of ADAFs and slim discs with a simple toy model. The advection 
‘cooling’ (per unit surface) term in the energy equation can be written as 


nadv 

4 2ttR? sU 


(see Eq. 239). 

Using the (non-relativistic) hydrostatic equilibrium equation 

R V K 

one can write the advection term as 


adv 


Q 


whereas the viscous heating term can be written as 


_ T «es c /m\ ( H_ 


2 R \rj J \R 


3 Kps C / TO 


Q+ = T--=- - , 


8 R \t] 


(108) 


(109) 


( 110 ) 
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where 


r=(^V 

\ K es R J 


Since £ a ~ 1, 


Q 


adv 




(112) 


(113) 


and, as said before, for geometrically thin discs ( H/R <C 1) the advective term Q adv is negligible compared to the 
heating term Q + and in thermal equilibrium viscous heating must be compensated by radiative cooling. Things are 
different at, very high temperatures, when (H/R) ~ 1. Then the advection term is comparable to the viscous term 
and cannot be neglected in the equation of thermal equilibrium. In some cases this term is larger than the radiative 
cooling term Q~ and (most of) the heat released by viscosity is a dvected toward the accreting body instead of being 
locally radiated away as happens in geometrically thin discs. 

From Eq.(57) one can obtain a useful expression for the square of the relative disc height (or aspect ratio): 


H\ 2 _V2(rh 


R 


K es V V 


(a£)' 


Using this equation one can write for the advective cooling 


Q adv = TUp a (aSp ( - 


The thermal equilibrium (energy) equation is 


Q+ = Q adv + Q~. 


(ff- 

(114) 

/n K . 


/ . \ 2 


TO \ 

kV/ 

(115) 


(116) 


The form of the radiative cooling term depends on the state of the accreting matter, i.e. on its temperature, density 
and chemical composition. Let us consider two cases of accretion flows: 

— optically thick 

and 

— optically thin. 

For the optically thick case we will use the diffusion approximation formula 


Q- = 8 -^ 
2 3k r E’ 


and assume «r = K es . With the help of Eq. (|114|) this can be brought to the form 

<■ • \ 1/2 

m \ 


Phick = 8T 


^es Rs 
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RsJ 




- 1/2 


V . 


For the optical thin case of bremsstrahlung radiation we have 

Q~ = 1.24 x 10 21 Hp 2 T 1/2 


which using Eq. (114) can be written as 

Qthin = 3-4 x 10~ 6 X 

• In the optically thick case we have therefore 
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FIG. 10: (a) Thermal equilibria for optically thick (The right solid 5-shaped line) and optically thin (the left solid line) accretion 
flows. The upper branches represent advection-dominated solution (ADAFs). Flows above the dotted lines r = 1 are optically 
thin - r is the effective optical depth calculated for radiation-pressure dominated (upper line) or gas-dominated (lower line) 
configurations. It is assumed that Mbh = IOMq, R = 5Rs, a = 0.1 and £ a = 1. (b) The same for a = 0.01.[From [4]] 


• In the optically thin case the energy equation has the form 

2 / t~i \ 1/2 


M- 
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0.18 


—X 
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, / m . 

aS vn + 
R V 


+3 x 10 " b a~ 2 ( — J (a£f = 0 


( 122 ) 


There are two distinct types of advection dominated accretion flows: optically thin and optically thick. We will 
first deal with optically thin flows known as ADAFs. 


1. Optically thin flows: ADAFs 


For prescribed values a and £ a , Eq. (122) is a quadratic equation in (■ m/rj ) whose solutions in the form of rh(E) 
describe thermal equilibria at a given value of R. Obviously, for a given £ this equation has at most two solutions. 
The solutions form two branches on the m(a£) - plane: 


the ADAF branch 


fn = 0.53k 6S v ( ) 

'\RsJ 


1/2 


C 1 ^. 


(123) 


and 


3/2 


• the radiatively-cooled branch 


fn = 1.9 x 10 5 r] ( —— 
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(124) 
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(1221 exists. This imp’ 


From Eqs. (1231 and (1241 it is clear that there exists a maximum accretion rate for which only one solution of Eq. 


ies the existence of a maximum accretion rate at 


m n 


1.7 x 10 3 ?7 a 2 


-1 

Rs) 


1/2 


(125) 


This is where the two branches formed by thermal equilibrium solutions on the m(aS) - plane meet as seen on Figure 

EH 

The value of m max depends on the cooling mechanism in the accretion flow and free-free cooling is not a realistic 
description of the emission in the vicinity (( R/Rs ) < 10 3 ) of a black hole. The flow there most probably forms a 
two-temperature plasma. In such a case m max ~ 10a 2 with almost no dependence on radius. For larger radii rii max 
decreases with radius. 


2. Optically thick flows: slim discs 


Since the first two terms in Eq. (1211 are the same as in (Eq. 122), the high to, advection dominated solution is 
the same as in the optically thin case but now represents the 


Slim disc branch 


m = 0.53 


B \ 1 / 2 

is) c ' oE - 


(123) 


Now, the full equation (121) is a cubic equation in to 1 / 2 and on the rh(aE) plane its solution forms the two 
upper branches of the S-curve shown in Fig. |10[ The uppermost branch corresponds to slim discs while the 
branch with negative slope represents the Shakura-Sunayev solution in the regime a. (see Sect. Ill F), i.e. 


• a radiatively cooled, radiation-pressure dominated accretion disc 


to = 160 Kgg 1 rj 



(aS)" 1 


(126) 


3. Thermal instability of radiation-pressure dominated discs 


Radiation-pressure dominated (P = P rac i) accretion discs are thermally unstable when opacity is due to electron 
scattering on electrons. Indeed 


d lnT 4 f 


eff 


din T r 


= 4 


because kr = K e s = const., while in a radiation pressure dominated disc Q + ~ vYh ~ HT 4 ~ T 8 /E so 
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(127) 


(128) 


and the disc is thermally unstable. This solution is represented by the middle branch with negative slope (see Eq. 
126) in Fig. 10 The presence of this instability in the model is one of the unsolved problems of the accretion disc 


theory because it contradicts observations which do not show any unstable behaviour in the range of luminosities 
where discs should be in the radiative pressure and electron-scattering opacity domination regime. 


4- Slim discs and super-Eddington accretion 


From Eqs. (114) and (126) one obtains for the disc aspect ratio 


H 
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= 0.11 


Rs 
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(129) 
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which shows that the height of a radiation dominated disc is constant with radius and proportional to the accretion 
rates. 


But this means that with increasing m advection becomes more and more important (see e.g. Eq. 113) and for 

(130) 


dr „ „ R 
- « 9.2 — 
V Rs 


advection will take over radiation as the dominant cooling mechanism and the solution will represent a slim disc. 
Equation (130) can be also interpreted as giving the transition radius between radiatively and advectively cooled disc 


for a given accretion rate m: 


Rt 


0.1 


Rs 


(131) 


Another radius of interest is the trapping radius at which the photon diffusion (escape) time Hr/c is equal to the 
viscous infall time R/v r 


i?, 


Hrv r HnYi M 


trapp 


= ^(-)R s . 
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(132) 


Notice that both i? trans and 7? trapp are proportional to the accretion rate. 

In an advection dominated disc the aspect ration H/R is independent of the accretion rate: 


H 

X =0.86£ a 


-) 
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(133) 


therefore contrary to radiatively cooled discs, slim disc do not puff up with increasing accretion rate. 
Putting (133) into Eq. (132) one obtains 


^=0.86C' ,! (|d 


1/4 


(134) 


Radiation inside the trapping radius is unable to stop accretion and since i?trapp ~ rn there is no limit on the accretion 
rate onto a black hole. 


The luminosity of the toy-model slim disc can be calculated from Eqs. (118) and (123) giving 

Q~ = oT 3 4 ff = 


n4 0.1 IvEdd 


& R 2 ’ 

which implies T e g ~ l/R}! 2 . The luminosity of the slim-disc part of the accretion flow is then 
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where we used Eq. (131). 


Therefore the total disk luminosity 

Etotal = -Idhin T l^slim = 


(135) 
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where Lthin is the luminosity of the radiation-cooled disc for which Eq. (581 applies. 


It is easy to see that the same luminosity formula L « LEdd(l + lnm) is obtained when one assumes mass-loss from 
the disc resulting in a variable (with radius) accretion rate: M ~ R. 

At very high accretion rates the disc emission will be also strongly beamed by the flow geometry so that and 
observer situated in the beam of the emitting system will infer a luminosity 


L sp h = -^Le dd(l + lnm), 


(138) 


where b is the beaming factor (see m for a derivation of b in the case of Ultra-Luminous X-ray sources). 

Numerical simulations do not seem to correspond to this analytical solutions (see e.g. Isa, m and [52] but they 
also disagree between themselves. The reasons for these contradictions are worth investigating. 

Additional reading: References Q], 0j, g], [SSL EE], @S], [5D], and [HD] ■ 
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VI. ACCRETION DISCS IN KERR SPACETIME. 

In this section we will present and discuss the set equations whose solutions represent a-accretion discs in the Kerr 
metric. This section is based on references 0 , Esa and 09] and to be understood requires some basic knowledge of 
Einstein’s General Relativity. 


A. Kerr black holes 


The components gij of the metric tensor with respect to the coordinates (t, x a ) are expressible in terms of the lapse 
N , the components p a of the shift vector and the components 7 a/ 3 of the spatial metric: 

dx l dx 2 = — N 2 dt 2 + ^ a p(dx a + P a dt)(dx& + P^dt), (139) 

which is a modern way of writing the metric. 

Remark 3. In this section only I will use conventions different from those used in other parts of the Chapter. 
First, I will use the so-called geometrical units that are linked to the physical units for length, time and mass by 


length in physical units 
time in physical units 

mass in physical units 


length in geometrical units, 

- length in geometrical units, 
c 

c 2 

— length in geometrical units. 


(140) 


Second, the radial coordinate will be called “r” and not “R". This should not confuse the reader since R is used 
only in the non-relativistic context where it denotes a radial coordinate and a radial distance, while in the relativistic 
context it only a coordinate. 


1. General structure , Boyer-Lindquist coordinates 


The Kerr metric in the Boyer-Lindquist (spherical) coordinates t, r, 9 : <p corresponds to: 
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y/AA’ 


p r = p* = 0, p* = -u, 
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with 


? = r 2 + a 2 cos 2 9 , A = r 2 — 2 Mr + a 2 , 
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where M is the mass and a = J/M is the angular momentum per unit mass. In applications one often uses the 
dimensionless “angular-momentum” parameter a* = a/M. 

Therefore in BL coordinates the Kerr metric takes the form of 


ds 2 


C 2 A , 9 A sin 2 9 . , , , 2 s 2 , 7 „ 

- dt 2 H--— (dip — ivdt) + —dr 2 + q 2 d9 2 . 


A 




A 


(145) 


The time (stationarity) and axial symmetries of the metric are expressed by two Killing vectors 


= e = v 


(v) ’ 


where is the Kronecker delta. 


(146) 
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Remark 4. Using Killing vectors (146) one can define some useful scalar functions: the angular velocity of the 
dragging of inertial frames w, the gravitational potential <f>, and the gyration radius 91, 


w = -rf e- 24, 9t 2 = — C~C- (147) 

In the Boyer-Lindquist coordinates the scalar products of the Killing vectors are simply given by the components of 
the metric, 


V ' V — 9tt, V ■ £ — <?iys> £ ' £ — 


(148) 


and therefore quantities defined in Eq. (147) can be explicitly written down in terms of the Boyer-Lindquist coordinates 


as: 


m 2 
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(149) 


• The horizon 

The black hole surface (event horizon) is at 


th = M + \J M 2 — a 2 . 


(150) 


Therefore a horizon exists for a* < 1 only. At the horizon the angular velocity of the dragging of inertial frame is 
equal to 


ujh = Hr 


a 

2Mr H ’ 


(151) 


where fin is the angular velocity of the horizon, i.e. the angular velocity of the horizon-forming light-rays with respect 
to infinity. The horizon rotates. 

The area of the horizon is given by 


S = 8nMrn = 87 tM \J M 2 — a 2 . 

The extreme (maximally rotating) black hole corresponds to 

a = M. 


(152) 


(153) 


For a > M the Kerr solution represents a naked singularity. Such singularities would be a great embarrassment 
not only because of their visibility but also because the solution of Einstein equation in which they appear violate 
causality by containing closed time-like lines. The conjecture that no naked singularity is formed through collapse of 
real bodies is called the cosmic censorship hypothesis (Roger Penrose). 

Remark 5. Rotation of astrophysical bodies 

Since this is a lecture in astrophysics let us leave for a moment the geometrical units. They are great for calculations 
but usually useless for comparing their results with observations. In the physical units 
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and therefore the maximum angular momentum of a black hole is 
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This is slightly more than the angular momentum of the Sun (Jq = 1.63 x 10 48 g cnr 2 s _1 , a® = 0.185): the gain in 
velocity is almost fully compensated by the loss in radius. 

For a millisecond pulsar which is a neutron star with a mass of ~ 1.4 Mq and radius ~ 10km the angular momentum 


is 


Jns = Ins^S 


8.6 x 10 48 
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(156) 


where Ins ~ o(x)M^sR^s 4S the moment of inertia and x = (MNs/Mg)(km/RNs) the compactness parameter. For 
the most compact neutron star x < 0.24 and a(x) < 0.489. Therefore for neutron stars that rotate at millisecond 
periods 


By definition 
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• the specific (per unit rest-mass) energy is 


£:=-rf-u, (158) 

• the specific (per unit rest-mass) angular momentum 

£:=£-u (159) 

and 


the specific (per unit mass-energy) angular momentum (also called geometrical specific angular momentum) 

3:=-§ = -i : l (160) 

Ct rj ■ u 


B. Privileged observers 

Let us consider observers privileged by the symmetries of the Kerr spacetime. The results below apply to any 
spacetime with the same symmetries, e.g. the spacetime of a stationary, rotating star. The four-velocity of a privileged 
observer is the linear combination of the two Killing vectors: 

u = z[r} + f2 obs £*) (161) 

where the redshift factor Z is (from the normalization u ■ u = 1) 

Z- 2 = ff- fj + 21W • S + fiobsl*• £ (162) 

Since for a ^ 0 the Kerr spacetime is stationary but not static, i.e; the timelike Killing vector 77 is not orthogonal 
to the space-like surfaces t =const. In such a spacetime ’’non-rotation” is not uniquely defined. 

Stationary observers are immobile with respect to infinity; their four-velocities are defined as 

<tat = (VV)~ 1/2 rf (163) 

but are locally rotating: £ sta t = & u st at 7 ^ 0 . 

The four-velocity of a locally non-rotating observer is a unit timelike vector orthogonal to the space-like surfaces 
t = const.: 

«zamo = e<S> Of + . (!64) 

defines four-velocity of the local inertial observer or ZAMO, i.e. Zero Angular Momentum Observers since 

•CzAMO = Ci u ZAMO = 6- 

Finally, in presence of matter forming a stationary and axisymmetric configuration, there are privileged observers 
comoving with matter. 
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TABLE I: Summary of properties of privileged observers 


Observer 

Four-velocity 

Angular velocity with respect 
to stationary observers 

Stationary 

u = (ff ■ rf)- 1 2 fj 

^stat = 0 

ZAMO (LNR) 

u = e"* (rj + ut) 

^ZAMO = W 

Comoving (with matter) 

u = Z (jj + 12£^ 

“com — 


C. The ergosphere 

For ZAMOs = ui but for stationary observers S7 — u = —ui. Therefore ZAMOs rotate with respect to infinity (but 
are locally non-rotating). They may exist down to the black hole horizon, where they become null: Uzamo -'Uzamo = 0. 

Stationary observers immobile with respect to infinity but rotating with angular velocity — u> with respect to ZAMOs 
can exist (their four-velocity must be timelike, 77-77 < 0 ) only outside the stationarity limit whose radius is defined 
by ff ■ ff= 0 : 

r ei (9) = M + \JM 2 - a 2 cos 2 9. (165) 

The stationarity limit is called the ergosphere. 


D. Equatorial plane 


We will discuss now orbits in the equatorial plane, where they have the axial symmetry. We are introducing the 
cylindrical vertical coordinate z = cosd is defined very close to the equatorial plane, 2 : = 0. The metric of the Kerr 
black hole in the equatorial plane, accurate up to the ( z/r)° is 


ds 2 = — 7 —^-dt 2 + ^ (dip — eodt) 2 + r 2 ^r 2 + dz 2 , 

A r 2 v ' d 

where now 

A = r 2 — 2Mr + a 2 , A = (r 2 + a 2 ) 2 - Aa 2 , co=^^, 

or simpler 

ds 2 = — (1 — dt 2 — 2 ujdtdp + ^dp 2 + -r-di' 2 + dz 2 . 

\ r J r z A 


(166) 


(167) 


(168) 


1. Orbits in the equatorial plane 

The four velocity of matter u l has components id, u v , u r , 

v? =u t 5\ t) +u*8 i w +u r 8\ r) . (169) 

The angular frequency with respect to a stationary observer, and the angular frequency 12 with respect to a local 
inertial observer are respectively defined by 

u v 

12 =—-, 12 = 12 — ui, 
ir 


( 170 ) 
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The angular frequencies of the corotating (+) and counterrotating (-) Keplerian orbits are 
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r 3 / 2 ± aM 1 / 2 ’ 


the specific energy is 
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r 2 — 2 Mr ± a(Mr) 1 / 2 
r (r 2 — 3 Mr ± 2o(Mr) 1 / 2 ) 1 ^ 


and the specific angular momentum is given by 

± ^ (Mr) 1 / 2 (r 2 =F 2o(Mr) 1 / 2 + a 2 ) 

r (r 2 — 3Mr ± 2a(Mr) 1 / 2 ) 1 ^ 2 
or 

(Mr) 1 / 2 (r 2 =F 2a(Mr) 1 / 2 + a 2 ) 
^ r 2 — 2Mr ± a(Mr) 1 / 2 


(171) 


(172) 


(173) 


(174) 


Both d and £ have a minimum at the last stable orbit, more often called ISCO (Innermost Stable Circular Orbit). 


Because of the rotation of space there is no direct relation between angular momentum and angular frequency but 


fj-t + ng-g m 2 Q-c c 
ff ■ f) + Qff ■ t; — 1 


(175) 


For the Sclrwarzschild solution (a = uo = 0) 


Zk = K 2 n K , (176) 

so a Newtonian-like relation (justifying the name “gyration radius” for 91) between angular frequency and angular 
momentum exists for (j. No such relation exists for £. 

ISCO 

The minimum of the Keplerian angular momentum corresponding to the innermost stable circular orbit (ISCO) is 
located at 


r ISCO 

Zi 

Z 2 


M {3 + Z 2 T [(3 - Z x )(3 + Z ± + 2Z 2 )] 1 ' 2 }, 


1 + (l - c 
(3a 2 /M 2 


2 /M 2 


. 1/3 


(1 + a/M) 1 / 3 + (1 — a/M) 1 / 3 


z t 


21 1/2 


Binding energy 

The binding energy 


at the ISCO is 

• 1 - vW9 « 0.06 for a = 0 

• 1 - /V3 « 0.42 for a = 1. 


(177) 


^bind — 1 /\ 


(178) 
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FIG. 11: Radii of characteristic orbits in the Kerr metric as a function of a* = a/M. The innermost stable circular orbit: 
r'isco, the marginally bound orbit riBCO (marked r m b), the photon orbit: r p h and the black hole horizon: vh (marked rh). 
(Courtesy of A. S§,dowski.) 


This corresponds to the efficiencies of accretion in a geometrically thin (quasi-Keplerian) disc around a black hole. 
For a Scliwarzschild black hole the frequency associated with the ISCO at nsco = 6M is 


Ffir(nsco) = 2197 



(179) 


IBCO 

The binding energy of a Keplerian orbit 1 — (£ = 0 at the marginally bound orbit (or IBCO: Innermost Bound 
Circular Orbit) 


„± 


= 2M T a + 2 \/M 2 =f aM. 


1 IBCO 

For a non-rotating black-hole rieco = 4M and the frequency associated with the IBCO is 

-l 


vk(ti bco) = 4037 ( — 


M 

■0 


Hz. 


ICO (Circular photon orbit) 

The Innermost Circular Orbit (ICO), i.e. the circular photon orbit is at 


r ph = 2 M ^ 1 + cos 

For a non-rotating black hole r p h = 3M. 


- cos 
3 


(w) 


(180) 


(181) 


(182) 


2. Epicyclic frequencies 

We will consider now consider a perturbed orbital motion in, and slightly off the equatorial plane. In the Newtonian 
case the angular frequency of such motions must be equal to the Keplerian frequency fix since in there is only one 
characteristic scale defined by the gravitational constant G. In General Relativity the presence of two constants G 
and c imply that the epicyclic frequency does not have to be equal to f Ik- 
The four-velocity for the perturbed circular motion can be written as 

u l = (l,u r ,u e ,n K + u v ) , (183) 

where u a are the velocity perturbations. 












34 


• For perturbations in the equatorial plane the equation of motion is 

d 2 


dt 2 


= 0, 


where 


= n 


r 2 — 6Mr ± 8 aM 1 / 2 r 1 / 2 — 3a 2 


K " 


(184) 


(185) 


. * 2 = 

at ISCO. In the Newtonian limit the epicyclic frequency equal to the Keplerian frequency k = f Ir¬ 


is the (equatorial) epicyclic frequency. In the Schwarzschild case a = 0 this is k 2 = 0^.(1 — 6 M/r) and vanishes 


• For vertical perturbations the equation is 


O r\ \ 

m +nK d^j ^ = ~ n±66, 

where the vertical epicyclic (angular) frequency is given by 

2 2 r 2 — 4aM 1 / 2 r 1//2 — 3 a 2 M 2 

n l = n K - -2 - 


(186) 


(187) 


In the Schwarzschild case (a = 0) the vertical epicyclic frequency is equal to the Keplerian angular frequency 
fix-, which is to be expected from the spherical symmetry of this solution. 

The angular velocity f2j_ appears also in the equation of vertical equilibrium of a (quasi)Keplerian disc which 
will be discussed later (sect. VIIIE|. Here let us just notice that Eq.(222) can be written as 


dp 2$o2 

-=-pe 


All these characteristic frequencies can be put into the form 


ft = f(x,a *) —, 


where x = r/M. For all relativistic frequencies x = a;(a*) and therefore they can be written as 


^ = F ( a *^ JJ- 


(188) 


(189) 


(190) 


Additional reading: Reference [2]. 


VII. ACCRETION FLOWS IN THE KERR SPACETIME 
A. Kinematic relations 

In the reference frame of the local inertial (non-rotating) observer the four velocity takes the form, 

U* = 7 (u l ZAMO + V^T 1 ^ + 


(191) 


The vectors t - '^ and are the unit vectors in the coordinate directions ip and r. The Lorentz gamma factor 7 
equals, 


7 = 


1 — (vivl) — (uM)^ 


(192) 
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The relation between the Boyer-Lindquist and the physical velocity component in the azimuthal direction is, 

= m, (193) 

which justifies the name of R - gyration radius. It is convenient to use the (rescaled) radial velocity component V 
defined by the formula, 


V 

Vi-v 2 

The Lorentz gamma factor may then be written as, 


= r yv (r) = u r g x J 2 . 


7 2 = 


1 


)(r 


1 


i -n 2 R 2 J V 1 ^ 2 


(194) 


(195) 


which allows writing a simple expression for V in terms of the velocity components measured in the frame of the local 
inertial observer, 


V = 


v ( r ) 

\J 1 — (u^)) 2 


v (r) 

\J\ - R 2 n 2 ' 


(196) 


Thus, V is the radial velocity of the fluid as measured by an observer corotating with the fluid at fixed r. 

Although a different quantity could have been chosen as the definition of the “radial velocity”, only V has directly 
three very convenient properties, all guaranteed by its definition: 

• (i) everywhere in the flow \V\ < 1, 

• (ii) on the horizon \V\ = 1, 

• (iii) at the sonic point |Vj « c s , 

where c s is the local sound speed. 

To see that property (i) holds, let us define 

V 2 = u r u r = u r u r g rr > 0. (197) 

Then, one has 

V 2 = V 2 /(l + V 2 ) < 1. (198) 

Writing V = \J~r 2 u r u r /(r 2 u r u r + A) demonstrates property (ii) since \V\ = 1 independent of the value of r 2 u r u r . 
For the proof of property (iii) of V see [3] . 

Other possible choices of the “radial velocity” such as u = |M r | are not that convenient. 


B. Description of accreting matter 

The stress-energy tensor T lk of the matter in the disk is given by 

T ik = (e+p) u l u k +pg ik + S lk + u k q l + u l q k , (199) 

where e is the total energy density, p is the pressure, 

Sik = vp<Ji k , ( 200 ) 

is the viscous stress tensor, p is the rest mass density and q l is the radiative energy flux. In the last equation v is the 
kinematic viscosity coefficient and is the shear tensor of the velocity field. From the first law of thermodynamic 
it follows that 


de = - dp + pT dS , 


P 


( 201 ) 
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where T is the temperature and S is the entropy per unit mass. Note, that in the physical units e = pc 2 + II, where 
II is the internal energy. For non-relativistic fluids, II <C pc 2 and p <C pc 2 , and therefore 

e + p~c 2 p. ( 202 ) 

We shall use this approximation (in geometrical units e + p ss p) in all our calculations. This approximation does 
not automatically ensure that the sound speed is below c, and one should check this a posteriori when models are 
constructed. We write the first law of thermodynamics in the form: 

dU = -pd + TdS, (203) 

where U = 11 /p. 


VIII. SLIM DISC EQUATIONS IN KERR GEOMETRY 

General-Relativistic effects play an important role in the physics of thin (H/r <C 1) accretion discs close to the 
black hole but they determine the properties of slim (H/r < 1) discs. We will derive the slim-disc equation and before 
discussing their properties we will say few words about thin discs. 

It is convenient to write the final form of all the slim disk equations at the equatorial plane, z = 0. Only these 
equations which do not refer to the vertical structure could be derived directly from the quantities at the equatorial 
plane with no further approximations. All other equations are approximated — either by expansion in terms of the 
relative disk thickness H/r, or by vertical averaging. 


A. Mass conservation equation 


From general equation of mass conservation, 


V 4 (p Ui ) = 0, 


and definition of the surface density E, 

r+H(r) 

E = / p(r, z)dz ~ 2Hp, 
J-H(r ) 

we derive the mass conservation equation, 

M = -2ttA 1/2 E — V 

In the Newtonian limit the mass conservation equation is: 


M = — 2-7T S« r . 


(204) 


(205) 


(206) 


(207) 


B. Equation of angular momentum conservation 

From the general form of the angular momentum conservation, 

v fc (T fci eo = o, 

we derive, after some algebra, 

2nr dr r dr \ r 4 dr) 

where F~ = 2 q z is the vertical flux of radiation, and 



(208) 


(209) 


£ = — (u£) = -u v = 7 


7*3 Al/2 


( 210 ) 
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is the specific (per unit mass) angular momentum. The term F~C represents angular momentum losses through 
radiation. Although it was always fully recognized that angular momentum may be lost this way, it has been argued 
that this term must be v ery small. Rejection of this term enormously simplifies numerical calculations, because with 
F~ £ = 0 equation (209) can be trivially integrated, 


-*>) = - 


r 4 dr 27t 


( 211 ) 


where £ 0 is the specific angular momentum of matter at the horizon (A = 0). In the numerical scheme for integrating 
the slim Kerr equations (with F~£ assumed to be zero) the quantity £o plays a n im portant role: it is the eigenvalue 
of the solutions that passes regularly through the sonic point. The rhs of Eq. (211) 1 represent the viscous torque 
transporting angular momentum. 

In the Newtonian case, a geometrically thin disc is Keplerian f2 m f \ K (see Eq. (214), £,x = £k = 7? 2 S2a', Eq. (211) 
takes the familiar form of 


i/E = 


M 

3tt 


1 - 


Ik 


(see Eq|57|. 


C. Equation of momentum conservation 

From the r-component of the equation V iT lk = 0 one derives 

V dV __ A 1 dP 

1 — V 2 dr r E dr ’ 

where P = 2 Hp is the vertically integrated pressure and 

MA (fl — f2^)(n — Q k ) 


A = - 


r 3 Afl+Ofc 


1 - Q 2 R 2 


( 212 ) 


(213) 


The Newtonian limit of Eq. (212) is 


Note that in Eq. (|212|) the viscous term has been neglected. 

Vr ^r-( Q 2 - Q K) 


r + — = 0 
r 


(214) 


For a thin disc: H/r ~ c 2 s /r£lx ^ 1 Eq. (214) is simply SI ss i.e. a thin Newtonian disc is Keplerian. The 
thickness of the disc depends on the efficiency of radiative processes: efficient radiative cooling implies a low speed of 
sound. 


D. Equation of energy conservation 

From the general form of the energy conservation 


v, (: T ik Vk ) = o, 

and the first law of thermodynamics, 





the energy equation can be written in general as 

Q adv 


Q + -Q~, 


where 





2 


(215) 


(216) 


(217) 


( 218 ) 
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is the surface viscous heat generation rate, Q is the radiative cooling flux (both surfaces) which is discussed in 
Section 5, and Q adv is the advective cooling rate due to the radial motion of the gas. It is expressed as 


T.V IS. 1 / 2 ^dS = M ^dS 
sj l — V 2 t dr 2nr dr 


(219) 


In stationary accretion flows advection is important only in the inner regions close to the compact accretor. In the 
rest of the flow the energy equation is just 


Q + = Q~. 


In the newtonian limit of Eqs. (218) and (211) one obtains 

3 GMM 


Q + = 


8tt R 3 


1 - 


( 220 ) 


(58) 


E. Equation of vertical balance of forces 


The equation of vertical balance is obtained by projecting the conservation equation onto the 9 direction 

hg VkTi = 0, where h l g = Sg — u x ug (221) 

and neglecting the terms 0 3 (cos9). For a non-relativistic fluid this leads to 

dP £ 2 - « 2 (£ 2 - 1) 

— = -pg z z = -p - - 4 -2 • (222) 


In the newtonian limit Eq. (222) becomes 


dP _ Ik 
dz ^ r 4 


(see Eq. 18). 


IX. THE SONIC POINT AND THE BOUNDARY CONDITIONS 


A. The “no-torque condition” 


There have been a lot of discussion about the inner boundary condition in an accretion disc. The usual reasoning 
is that for a thin disc the inner boundary is at ISCO and since it is where circular orbits end the boundary condition 
should be simply that the “viscous” torque vanishes (there is no orbit below the ISCO to interact with). Several 
authors have challenged this conclusion but a very simple argument by Bohdan Paczyriski [33] shows the fallacy of 
these challenges. 


Using Eq. (206), one obtains from Eq. (211) 


A 3 / 2 7 2 1 dU 

4 £ — £q dr 


(223) 


Next, from the viscosity prescription v « and taking for simplicity the non-relativistic approximation (this 

does not affect the validity of the argument but allows skipping irrelevant in this context multiplicative factors) one 
can write 


a H z 


dil 


' — in dr 


— kuH z 


£-4 r 


n 


(224) 


where v v = Rfl Although we have dropped the GR terms, the equation (2241 does not assume that the radial velocity 
is small, i.e.this equation holds within the disk as well as within the stream below the ISCO. 
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Far out in the disk, where l to, one obtains the standard formula (see Eq. 36) 


V r 


a v u 


R > Ri- 


(225) 


The flow crosses the black hole surface at the speed of light and since it is subsonic in the disc it must somewhere 
become transonic, i.e. to go through a sonic point, close to disc’s in ner ed ge. 

At the sonic point we have v r = c s ~ (H/R)v v , and the equation (224) becomes: 


= 1 


6 . 

•°in 




R\n f’in — A) 


R = R\, 


(226) 


If the disk is thin, i.e. H ln /R ln <C 1, and the viscosity is small, i.e. «<1, then Eq. (226) implies that (4-4)An -C 1, 
i.e. the specific angular momentum at the sonic point is almost equal to the asymptotic angular momentum at the 
horizon. 

In a steady state disk the torque X has to satisfy the equation of angular momentum conservation (211), which can 
be written as 


X = M (t - 4), X m = M {t in - to) ■ (227) 

Thus it is clear that for a thin, low viscosity disk the ‘no torque inner boundary condition’ (X ln « 0) is an excellent 
approximation following from angular momentum conservation. 

However, if the disk and the stream are thick, i.e. H/r ~ 1, and the viscosity is high, i.e. a ~ 1, then the angular 
momentum varies also in the stream in accordance with the simple reasoning presented above. However, the no-stress 
condition at the disc inner edge might be not satisfied. 

Additional reading: Reference [B], 
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Appendix 


Thermodynamical relations 


The equation of state can be expressed in the form: 


71 11 B 2 

P = P r + -pTi + — pT e + —, 

Pi Pe 247T 


(228) 


where P r is the radiation pressure, 7 Z is the gas constant, p; and p e the mean molecular weights of ions and elec¬ 
trons respectively, T), and T e ion and electron temperatures, a the radiation constant (not to be confused with the 
dimensionless angular momentum a in the Kerr metric), and B the intensity of a isotropically tangled magnetic field, 
includes the radiation, gas and magnetic pressures. The radiation pressure P r , the gas pressure P g , and the ma gneti c 
pressure P m correspond respectively to the first term, the second and third terms, and the last term in equation (2281. 

The mean molecular weights of ions and electrons can be well approximated by: 


Pi 


4X + Y 


l + X’ 


(229) 
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where X is the relative mass abundance of hydrogen and Y that of helium. We may define a temperature as 


T = n ("— + — ) , 

\Pi Me/ 


where 


, 1 1 
M —-1- 


Hi p e J 1 + 3X + 1/2 Y 


(230) 


(231) 


is the mean molecular weight. In the case of a one-temperature gas (T) = T e ), one has T = T) = T e . For an optically 
thick gas, P r = (4cr/3c)T/. 

For the frozen-in magnetic field pressure P m ~ P 2 ~ p 4//3 , therefore we may write the internal energy as 


U = — Tt 


nr 


e a p 


1/3 


pc r pm u {'y g - 1) 

where e 0 is a constant (P m = l/3e oy o 4 / 3 ) and 7g is the ratio of the specific heats of the gas. We define 


P P 4 — ft 

P = ^, ^ = = 

P Pg H” Pm 3p m 


(232) 


(233) 


From equations (2281 and (232) one obtains the following formulae (see e.g. “Cox & Giuli” 2004) for the specific heat 
at constant volume: 


c v = 


n 


M(7 g - 1) 


12(1 — P/Pm){lg — 1) + /3 

P 


4 — 3/3* P 

r 3 -i ~pT 


and the adiabatic indices: 


F, - 1 = 


(4-3/?*)( 7g -l) 

12(1 — P/Pm)(lg — 1) + /3 


(234) 


(235) 


r 1 =/3* + (4-3/3*)(r 3 -l). (236) 

The ratio of specific heats is 7 = c p /cv = T'i/P- For (3 = f3 m we have T 3 = 7g and F 3 = (4 — /3)/3 + /I( 7g — 1). For 
an equipartition magnetic field (/3 = 0.5) one gets r 3 = 1.5 and for /3 = 0.95, F-| = 1.65 (here we have used 7g = 5/3). 
One expects /3 m ~ 0.5 — 1. Since 


dS d In T / d In E dlnH\ 

T dR = Cv [~dR~ {Ts ~ 1} {-^r — dlTJ\ ’ 


The advective flux is written in the form: 


Q 


adv 


M P ^ 

mr-p 


where 


& 


4-3/?* dlnT 
r 3 - 1 dlnR 


+ (4-3/3*) 


dlnS 

dlnP 


(237) 


(238) 


(239) 


The term oc din H/d In R has been neglected. Since no rigorous vertical averaging procedure exists, the presence 
or not of the d In H/dhi R - type terms in this (and other) equation may be decided only by comparison with 2D 
calculations. 

The formulae derived in this section are valid for the optically thin case r = 0 if one assumes /3 = /3 m . 


Reference: Weiss, A., Hillebrandt, W., Thomas, H.-C., & Ritter, H. 2004, Cox and Giuli’s Principles of Stellar 
Structure, Cambridge, UK: Princeton Publishing Associates Ltd, 2004. 
























